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Whether a quantum bath can be approximated as classical noise is a fundamental issue in central 
spin decoherence and also of practical importance in designing noise-resilient quantum control. Spin 
qubits based on bismuth donors in silicon have tunable interactions with nuclear spin baths and 
are first-order insensitive to magnetic noise at so-called clock-transitions (CTs). This system is 
therefore ideal for studying the quantum/classical nature of nuclear spin baths since the qubit-bath 
interaction strength determines the back-action on the baths and hence the adequacy of a classical 
noise model. We develop a Gaussian noise model with noise correlations determined by quantum 
calculations and compare the classical noise approximation to the full quantum bath theory. We 
experimentally test our model through dynamical decoupling sequence of up to 128 pulses, finding 
good agreement with simulations and measuring electron spin coherence times approaching one 
second — notably using natural silicon. Our theoretical and experimental study demonstrates that 
the noise from a nuclear spin bath is analogous to classical Gaussian noise if the back-action of the 
qubit on the bath is small compared to the internal bath dynamics, as is the case close to CTs. 

However, far from the CTs, the back-action of the central spin on the bath is such that the quantum 
model is required to accurately model spin decoherence. 

PACS numbers: 03.65.Yz, 76.30.Mi, 76.60.Lz 


Introduction. Central spin decoherence due to coupling 
to the environment is not only a central issue in under¬ 
standing quantum-to-classical transitions mi, but also 
one of the key challenges in the realization of quantum 
computation [3]. There are two distinct models to de¬ 
scribe the decoherence processes in such cases: in the 
semiclassical model, the central spin accumulates ran¬ 
dom phases due to thermal or quantum fluctuations of 
the environment l 1, while in the quantum model, 
the coupling between the central spin and the environ¬ 
ment produces entanglement and results in leakage of 
the which-way information from the central spin to the 
environment EH3- The fundamental difference between 
these two models lies in the fact that the classical noise is 
independent of the central spin state while the quantum 
noise is governed by the back-action from the central spin 

mm- 

The classical noise model of quantum baths is a useful 
approximation in designing noise-resilient quantum con¬ 
trol m, which would otherwise require a large amount of 
numerical simulations of many-body dynamics of quan¬ 
tum baths. Dynamical decoupling has been employed to 
extract the noise spectra of baths j!3Hl6j . which are in 
turn used to design optimal quantum control for protect¬ 
ing quantum coherence and quantum gates. The viability 


of such methods critically depends on whether the noise 
picture is valid or not. Therefore, examining the condi¬ 
tions for the validity of classical noise model is not only 
of fundamental interest but also highly desirable for ac¬ 
curate quantum control under realistic conditions. 

Bismuth donors in silicon (Si:Bi) have recently at¬ 
tracted much attention in spin-based quantum compu¬ 
tation due to a number of favourable properties Ha¬ 
ng. These include long electron spin coherence times 
of up to three seconds [5(11 observed for Bi donors (in 
isotopically-enriched silicon-28) tuned to so-called clock 
transitions (CTs) — also known as optimal working 
points (OWPs) [21] or zero first-order Zeeman (ZEFOZ) 
transitions) — whose frequency is insensitive, to first or¬ 
der, to magnetic field fluctuations. The coherence times 
of donor electron spins in natural silicon is typically lim¬ 
ited to a few hundred microseconds by the 5% naturally 
abundant 29 Si nuclear spins. At CTs, the effect of the 
29 Si on the electron spin coherence is strongly suppressed 
(though not completely removed), leading to coherence 
times of up to 100 ms. Previous studies have focused on 
quantum approaches to model electron spin decoherence 
from 29 Si nuclear spin baths [221 GHj, however, perform¬ 
ing such calculations near the CTs can be challenging 
due to the strongly correlated nuclear spin baths [22]. 
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In this Letter, we explore the applicability of a classi¬ 
cal Gaussian stochastic noise description [24[ [25] of the 
nuclear spin bath around Bi donors in natural silicon, 
especially near the Si:Bi CTs. Such classical Gaussian 
noise is fully characterized by the two-point correlation 
functions or the noise spectrum [25] . We demonstrate 
the validity of such a semiclassical model by comparisons 
with exact results from the quantum model and with ex¬ 
perimental measurements. 

System and Hamiltonian. For the Si:Bi system inter¬ 
acting with a 29 Si nuclear spin bath (/,; = 1/2 and natural 
abundance 4.7% throughout the host lattice), the system 
Hamiltonian is divided into three parts [22] [23]: 

H = H cs + Hi n t + Hbath, (1) 

with 

H CB =w e S z -w* i I5+A 0 S-I 0) ( 2 a) 

flint =S Z J2 A i ( 2b ) 

i 

Hbath = -w® X '■ » % • h , (2c) 

i i<j 

where S is the Bi donor electron spin operator, I 0 (I*) is 
the 209 Bi ( 29 Si) nuclear spin operator, uj e , w ® 1 and w ® 1 
are correspondingly the Larmor frequencies of the donor 
electron spin, 209 Bi and 29 Si nuclear spins (which are 
related to their gyromagnetic ratios 7 by uj a = r ) a B) 1 
A 0 (Ai) is the coupling strength between the donor elec¬ 
tron spin and the 209 Bi( 29 Si) nuclear spins, By is the 
nuclear-nuclear interaction tensor, and B is the mag¬ 
netic field applied along the z axis. Here we have ne¬ 
glected the non-secular terms in H lnt which induce the 
central spin relaxation, because the qubit energy splitting 
is much larger than the qubit-bath coupling (see online 
Supplementary Materials for detailed discussions on the 
effects of the non-secular terms on the pure-dephasing 
model) [26]. 

Semiclassical model for quantum decoherence. The 
combined Bi donor electron-nuclear spin system (S = 
1/2, and Iq = 9/2) with the Hamiltonian H cs has 20 
eigenstates with eigenenergies dependent on the magnetic 
field [26]. By focusing on a pair of the eigenstates, |+) 
and |—), of the central spin, we can recast the system 
Hamiltonian in Eq. 0 as a function of these two central 
spin states: 

H (±) = ± P+ + P+ \ P ~ P Z + H bath , (3) 

where (3 Z = Pt = A-iPf is the Overhauser field 
operator and P± = (±|£ z |±). The CTs are charac¬ 
terised by P+ ~ P- such that central spin decoherence 
at the |+)•<-> | —) transition is strongly suppressed due to 
the nearly identical evolutions of the nuclear environment 
conditioned on the central spin state (i.e. H ~ H^'). 



FIG. 1: (Color online) (a) Relative two-point correlation func¬ 
tions C(t) — (7(0) at the CT (Bct = 79.9 mT) calculated 
by different orders of CCE. Here we choose a specific nu¬ 
clear spin configuration with B || [110]. (b) C(t) — (7(0) 

(solid lines) at the CT for several magnetic field orienta¬ 
tions in the [001] — [HO] plane with 0 = 0° corresponding 
to [001]. Results are obtained by averaging over 50 differ¬ 
ent nuclear spin configurations. Dashed lines are fits of the 
form A 2 {exp[—(Itl/r) 71 ] — 1}. (c) The fitting parameters as 
functions of 9. 


Consequently, the back-action of the central spin on 
the environment uni m is quite small near the CTs 
(|P + — P_| < |P + + P_|), so we may infer that a semi¬ 
classical model for quantum decoherence should well re¬ 
produce the results from the quantum model. 

The two-point correlation function of the nuclear spin 
noise is defined as 

C(t) = {/3 z (t)/3 z ( 0)) = (e iHet l3 z (0)e~ iHet p z (0)), (4) 

where (• • •) = Tr[p(, • ■ ■ ] denotes the ensemble average 
over the density matrix pb = E/2 M for M nuclear spins 
at infinite temperature, and H e = ^ P+ ^ P ~^ /3 Z + H hath is 
the effective Hamiltonian for the nuclear spin bath. We 
choose the factor 1 +l 2 ' instead of + 9 in H e for 
the following considerations: (i) near the CTs, P+ ~ P_, 
so these two factors differ only in their sign which has no 
effect on the nuclear spin dynamics; (ii) far away from 
the CTs, P+ S3 —P_ S3 1/2, so the chosen factor can 
better represent the back-action of the central spin on 
the nuclear spin bath. 

In dynamical decoupling (DD) control, a se¬ 
quence of N 7 r-flips are applied to the central 
spin at times {ti ,£2 • • • t/v} to suppress magnetic 
noise [271 ES]. In this paper we consider IV-pulse 
Carr—Purcell—Meiboom—Gill (CPMG-7V) control (2jT, 
150} with tk = (2 k — 1)/2IV (k = 1, 2, • ■ ■ ,N). With the 
Gaussian approximation [25] , the central spin decoher- 
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ence under DD control is 


L(t) = exp 


= exp 



f f dtidt-2C(ti 

Jo Jo 







(5) 


where P e = |P + — P_| is the scaling factor, and £7(w) = 
f_ e lult C(t)dt is the noise spectrum of the nuclear spin 
noise [since C(t) = £7(—t), so £7(w) = £7(—w)], f(t) = 
(— l) k for [tk,tk+i] is the modulation function, P(wt) = 
I Efc = o(- 1 ) fc ( c “ rtfc+1 — e tultk )\ 2 is the filter function. 

Characterization of the nuclear spin noise. The corre¬ 
lation functions C(t) for the natural 29 Si nuclear spin 
bath in silicon are calculated by the CCE method 
0 m ng. As a typical example, we choose to study 
the central spin transition between |+) = |5, —1) and 
|—) = |4, —2), using the basis |F, mj?) where F(= I + S) 
is the total spin and 77 Jf(= ms + mi) is its projec¬ 
tion. This transition is a CT at B = 79.9 rnT (where 
P + = P_= 0.0525). 

In Fig.^a), we show the relative correlation function 
C{t) — £7(0) corresponding to the quantum fluctuations 
of nuclear spin noise calculated for a random nuclear 
spin configuration with B || [110]. Since C(t) depends 
on the specific nuclear spin configuration, we show in 
Fig- 0b) the relative correlation functions averaged over 
many bath configurations. The results can be well fitted 
by a stretched exponential decay: 

£7(i) = C(0) + A 2 {exp[—(|t|/r)"] — 1}, (6) 


where A is the correlation amplitude, r is the correla¬ 
tion time and n is the stretch factor. In Fig. [TJc), we 
show the fitting parameters as functions of the mag¬ 
netic held orientation. As the held direction varies from 
[001] —► [111] —>• [110], the correlation amplitude A hrst 
increases with 6 (the angle from [001]), reaches a max¬ 
imum at about 55° ([111]) and then slightly decreases. 
The stretch factor n has about the same trend as A apart 
from some oscillations due to systematic fitting errors, 
while the correlation time r has the opposite trend. The 
dependence of the correlation functions on the magnetic 
held orientation is due to the anisotropy of the dipolar 
interaction between nuclear spins and the silicon lattice 
structure m Ea, and can be well understood from a 
microscopic analysis of the nuclear spin dynamics. 

As the major contribution [seeQa)], the pairwise hip- 
hop processes in the nuclear spin bath lead to a correla¬ 
tion function as 

! 2 Z 2 D 2 1 

Z 2 \ p2 [ C0S ( W D'*) - !] > , (7) 
'7 ij ) 


where = 2y Z 2 + D 2 - is the noise frequency, Zij = 
l P+ l^l P ~l (Ai_Aj) is the energy cost of hip-hop processes, 


(a) B C t + 0.15 mT 


(b) B ct + 9 mT 
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(c) B C j + 35 mT 


(d) High-field limit (470 mT) 
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FIG. 2: (Color online) Comparisons of electron spin de- 
coherence obtained by the quantum model (solid lines), the 
semiclassical model (dashed lines)] and the experimental mea¬ 
surement (circles). Measurements were made using (a-c) the 
15, — l)-f->-|4, — 2) transition at various fields shifted from the 
CT at 79.9 mT, or (d) the 15, — 4)e->j4, —5) transition close to 
the high-field limit at 468.65 mT. (e) The spin decoherence 
time for T 2 e under CPMG control obtained from the models 
and experiment for various fields. Here N= 2, 4, 8, 16, 32, 64, 
128 corresponds to the DD control CPMG-2, XY-4, XY-8, 
XY-16, (XY-16)x2, (XY-16)x4, (XY-16)x8. In theoretical 
calculations, CPMG-1V is equivalent to XY-DD. The exper¬ 
imental data is corrected to remove effects of instantaneous 
diffusion and spin relaxation. The theoretical results are ob¬ 
tained by averaging over 100 nuclear spin conhgurations. 


Dij = 7 2 (3cos 2 dij — l)/4|Rjj| 3 is the dipolar interaction 
strength, R,j is the displacement between the *th and 
jth nuclear spins, and dij is the angle between R,j and 
B. For a given Z. l;) . the correlation amplitude and noise 
frequency increase with the dipolar interaction strength 
Dij. When B || [001], the nearest-neighbor nuclear spin 
pairs have zero dipolar interaction, so the flip-flop pro¬ 
cesses mainly occur between the 2rd and 3rd-nearest 
neighbors. This produces a minimum in A and n and a 
maximum in r, as a result of the relatively slow nuclear 
spin dynamics. When B || [111], the nearest-neighbor 
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FIG. 3: (Color online) Calculated electron spin decoherence 
under an exact quantum model (solid lines) and semiclassical 
model obtained from noise spectroscopy of the CPMG-100 
DD (dashed lines), evaluated (a) far from the CT (Bct + 
1000 G) and (c) close to the CT (Bct + 10 G). (b) Filter 
function F(uot)/(u> 2 t) for CPMG-100 noise spectroscopy, with 
t = 10 ms. (d) Comparison of the noise spectrum from the 
CPMG-100 spectral decomposition in (c) to the exact one 
from CCE calculations. Here we choose a specific nuclear 
spin configuration with B || [110]. 


nuclear spin pairs have the strongest dipolar interaction, 
so A reaches a maximum, r a minimum, and n ss 1 since 
the relatively fast nuclear spin dynamics make the noise 
like a classical Lorentzian noise. 

Comparisons with experiments and quantum model. To 
explore the validity of the semiclassical model, we com¬ 
pare in Fig. [2] the decoherence obtained from Eq. ([ 5 ]) with 
both the results from the CCE method [8] and the ex¬ 
perimental measurements. The measurements were con¬ 
ducted on a nat Si sample doped with 209 Bi at a concen¬ 
tration of 3xl0 15 cm -3 . The experiments were realized 
at 4.8 K, which gives an electron spin relaxation time 
Ti e of 3.5 s. The magnetic field was aligned close to 
the [241] crystal axis. The DD 7 r pulses were applied 
through adiabatic fast passages (hyperbolic secant func¬ 
tions, 20 /jls in duration spanning 6 MHz close to the CT 
and 12 MHz in the high-field limit) in order to achieve 
high-fidelity operations despite the ESR linewidth, which 
ranged from 6 to 12 MHz at the fields studied here [25] , 
Owing to the high-fidelity DD control via the adiabatic 
passage method, remarkably, we managed to apply more 
than 100 control pulses and therefore extend the electron 
spin coherence times to one second in the natural silicon 
sample. We expect this can be further increased by using 
optimal magnetic field orientations and higher-order DD. 

The results from the quantum model and the experi¬ 
mental data are in good agreement for all the the mag¬ 
netic fields. However, this is not the case for the semi¬ 


classical model. Close to the CT [Hct + 0.15 mT in 
Fig. [2ja) , where P + = 0.0527, P_ = 0.0521], the semi¬ 
classical coherence, quantum coherence and experimental 
results coincide for various DD control sequences, indi¬ 
cating that the nuclear spin bath is well described by a 
classical noise. As the magnetic field is shifted away from 
the CT [Hot + 9 mT in Fig. [2](b), where P + = 0.0695, 
P- = 0.0340] the semiclassical model begins to show 
small deviations from the quantum model for DD with 
large numbers of pulses. Further still [Hot + 35 mT in 
Fig. [2](c), where P + = 0.1172, P_ = —0.0198], there 
are clear deviations between the semiclassical and quan¬ 
tum models for all levels of DD control, but the time 
scales remain nearly equal. Finally, for the transitions 
in the high-field limit [B = 470 mT in Fig. |2][d), where 
P + = 0.4264, P_ = —0.5], the semiclassical model ceases 
to be valid and shows significant differences from the 
quantum model and the experiment. These comparisons 
demonstrate the classical nature of the nuclear spin noise 
near the CT. This is understandable since the feedback 
of the central spin on the evolution of the nuclear spin 
bath (~ l P+ ~ P ~l ) is largely reduced at the CT. 

Limitations of DD noise spectroscopy method. Pre¬ 
vious studies have adopted the DD noise spectroscopy 
method to characterize the baths Bang. The main 
idea is to use a specific DD control sequence (such as 
CPMG-fV with large N) with the filter function approx¬ 
imated as a Dirac delta function at loq = ±nN/t [see 
Fig. |3](b)], i.e., F(ut)/(uj 2 t) « 1 r[6(w - u 0 ) + S(uj + w 0 )], 
and following Eq. (|5| to determine the noise spectra as 
S(±u 0 ) = — 21n[P(t)]/(fP 2 ). This method relies on the 
validity of the semiclassical model. 

However, the noise model may be insufficient when the 
back-action of the central spin on the environment dy¬ 
namics is significant. For example, in the nat Si:Bi system, 
the noise model is invalid for transitions far away from 
CTs. To demonstrate this point, we use the DD noise 
spectroscopy method to determine the effective noise 
spectra corresponding to the CPMG-100 case, and then 
use the derived noise spectra to calculate the spin deco¬ 
herence under other DD control sequences. In Fig. [3](a) 
we show the comparisons between the exact decoherence 
model and the semiclassical model using the DD noise 
spectroscopy method, and find increasing discrepancies 
as the pulse number of CPMG-fV deviates from 100. In 
contrast, close to the CTs, the DD noise spectroscopy 
method can not only reproduce the spin decoherence 
curves for other DD control [see Fig. |3][c)], but also well 
reproduce the exact noise spectrum obtained from CCE 
calculations [see Fig. |3jd)] . Here the decoherence profile 
shows increasingly violent oscillations (corresponding to 
the relatively high-frequency noise) for a single nuclear 
spin configuration, reducing the efficiency of DD control 
at long-time scales. The reason is that when the DD 
modulation frequency matches the noise frequency, the 
noise can be amplified rather than suppressed |33j . 
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Summary. We have presented a semiclassical model to 
study the decoherence of electron spin qubits in natural 
silicon near the CTs in nat Si:Bi system. The comparisons 
of the semiclassical results against the exact quantum 
results and experimental measurements demonstrate that 
the nuclear spin bath acts more and more like a classical 
Gaussian noise as the CTs are approached. Our findings 
deepen the understanding of spin baths near CTs and 
are useful for optimizing the DD control in silicon-based 
quantum computation indeed we have already shown 
here that using DD at CTs, electron spin coherence times 
of about one second can be measured in natural silicon. 
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